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Abstract 

Weak gravitational lensing from large-scale structure enhances and 
reduces the fluxes from extragalactic point sources with an r.m.s. am- 
plitude of order 15%. In cosmic microwave background (CMB) exper- 
iments, sources exceeding some flux threshold (/)c are removed, which 
means that lensing will modulate the brightness map of the remaining 
unresolved sources. Since this mean brightness is of order 100/xK at 
30 GHz for a reasonable flux cut, one might be concerned that this 
modulation could cause substantial problems for future CMB experi- 
ments. We present a detailed calculation of this effect and, fortunately, 
find that its power spectrum is always smaller than the normal point 
source power spectrum. Thus although this effect should be taken into 
account when analysing future high-precision CMB measurements, it 
will not substantially reduce the accuracy with which cosmological pa- 
rameters can be measured. 
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1 INTRODUCTION 



Future high-precision measurements of the fluctuations in the cosmic mi- 
crowave background (CMB) may allow key cosmological parameters (such 
as ri, A, ilfi, the Hubble constant, etc.) to be measured with unprecedented 
accuracy (Jungman et al. 1996; Bond et al. 1997; Zaldarriaga et al. 1997), 
and a number of ground-, balloon- and space-based missions are currently 
being planned for this purpose. Both when designing such missions and 
when analyzing the data sets that they produce, it is important that all rel- 
evant sources of foreground contamination are well understood, to minimize 
the risk that foreground signals are misinterpreted as CMB fluctuations. 
Since numerous experiments are are currently in the planning and design 
stages, it is therefore timely to catalog and quantify all possible foreground 
effects. Figure |l| summarizes recent foreground estimates from Tegmark 
& Efstathiou (1996), hereafter "TE96", and Bersanelh et al. (1996). The 
purpose of this Letter is to investigate yet another foreground effect which 
has not been previously studied: that of weak gravitational lensing of point 
sources. In other words, we will discuss the extent to which the point source 
region (which occupies mainly the lower right corner of Figure |l]) is expanded 
by lensing. 

Since it is impossible to remove all point sources (as there are for all 
practical purposes infinitely many, and not all with the same spectra), the 
standard procedure in CMB experiments is to remove all point sources below 
some flux cut (/)c, either by subtracting their estimated emission or by dis- 
carding all contaminated pixels. The total sky brightness from the remaining 
sources tends to be dominated by faint ones, whereas the brightness fluctua- 
tions (which contribute to Figure |l]) are dominated by the sources just below 
the flux cut. For this reason, the average brightness B exceeds the r.m.s. 
fluctuations by a substantial factor, typically ~ 10^, as shown in Figure 
A process that caused even minor spatial modulations of B could therefore 
pose a serious problem for CMB experiments. 

Weak gravitational lensing by large-scale structure (see Villumsen 1996 
for a recent discussion) can affect observed power spectra in more ways than 
one. If we imagine a pattern painted on the inside of a rubber (celestial) 
sphere, weak lensing corresponds to stretching and compressing the rubber 
in a random fashion, much like the way our mirror images get distorted by 
non-flat mirrors in amusement parks. For weak lensing, this distortion will 
always constitute a one-to-one mapping of the image, i.e., there will be no 
caustics where the image "folds over" on itself. There is thus no smearing 
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involved, so the fluctuation power is conserved. Instead, there is a "Robin 
Hood effect" where power is redistributed between multipoles, from those 
which have more to those which have less (Seljak 1996). For the CMB power 
spectrum, this cff'cct is typically a few percent, and may well be detectable 
in future CMB experiments. For point sources, on the other hand, this 
effect is completely absent: since they have a Poisson power spectrum with 
Q constant, the Robin Hood effect (which effectively simply smoothes the 
power spectrum), will leave Ci unaffected. 

However, the presence of a flux cut for point sources produce a different 
effect, which is absent for the CMB. In those regions where lensing causes 
magnification, we will see (and remove) disproportionately many sources 
above the flux cut (Turner et al. 1984, Broadhurst et al. 1995). Since lensing 
leaves the average brightness unaffected, the sky brightness from unremoved 
sources will be lower than average in this region. In other words, weak 
lensing causes the type of modulation of B that we warned about above. The 
rest of this Letter is organized as follows. We present a detailed calculation 
of this effect in Section 2, assess its importance in Section 3 and summarize 
our findings in Section 4. 



2.1 How magnification bias affects the source counts 

Let the unlensed source count function n(> 0) denote the number of sources 
per steradian whose flux would exceed (/) in the absence of lensing. When 
weak gravitational lensing magnifies a patch of sky with a magnification 
factor m, this has two separate effects on n(> 0): 

• The sources become a factor m brighter. 

• The sources are seen further apart, reducing their number density by 
a factor m. 

In summary, the lensed source count function, which we denote n(> 0), is 
given by 



2 CALCULATION OF THE EFFECT 




(1) 
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Defining the differential source count as = dn{< 4')/d(j) = —dn[> 

(t))/d(l) (a positive function), we tlius have 

^'(<^) = ii) ■ (2) 



For weak lensing, we can write m = 1 + A, where |A| <C 1 (Broadhurst et 
al. 1995), and we wili make the approximation of dropping aU quadratic and 
higher order terms in A. Using n'[(f)/m) ~ n'{(j) — (/)A) ~ n'{(j)) — n"{<f))(j)A, 
we can rewrite equation (^) as 

n'{(t>,r) = n'{<t>) - \ [n'(0)<A'l ' A(?), (3) 

where we have exphcitly indicated the fact that the magnification field de- 
pends on f , the unit vector pointing in our direction of observation. 

Since the total sky brightness /q°° n'cpdcj) is finite, the quantity n'cf)'^ must 
clearly approach zero both as (/> — > and as (/> — > oo, otherwise the integral 
would diverge at the faint or bright end. Since the second term in equa- 
tion (^) becomes a total differential when multiplied by (j), this means that 
the total brightness is not affected by lensing. This well-known fact is also 
verified by integrating (p times equation (|2|) and changing variables. How- 
ever, when the upper integration limit is not oo, the same procedure shows 
that the brightness contribution from all sources below some fixed flux cut 
is affected by lensing, and we will now evaluate this effect. 



2.2 The effect on the point source correlation function 

For CMB purposes, we can to a good approximation (TE96) model the 
locations of extragalactic point sources as completely uncorrelated. In the 
absence of lensing, we can therefore write the observed density of point 
sources above some flux cut (j)c as a sum of angular delta functions, n(r) = 
X]<^(?, fi)) ^-iid model n(r) as a Poisson process. A Poisson process satisfies 
(see e.g. Appendix A of Feldman, Kaiser & Peacock 1994) 

(n(f))p = rl(f), (4) 
(n(f)n(f'))p = n(f)n(f') -|- (5(f , r')n(f). (5) 

We use Op to denote ensemble averages with respect to the Poisson process, 
to distinguish these from ensemble averages with respect to the random field 
A, which we will denote by ( )/. When performing both averages below, we 
will omit subscripts and write ( ) = (( 
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Magnification bias forces us to take into account that the total source 
population is the union of a number of independent subpopulations, corre- 
sponding to different flux classes. We therefore write 

n(r) = / n{(j),r)d(j), (6) 



where n'{<j),r)d(j) is a sum of delta functions corresponding to the sources 

whose fluxes fall between (f) and (p + dep. With this notation, equations (§) 
and dll) become generalized to 

(n'(</.,?))p = fi'icP,?), (7) 

(n'(0,f)n'(0',?'))p = n'{^,^)n'{^',^') + 5{(P-^')6{^,?')n'{^,r). (8) 

The observed sky brightness a; in a direction f is clearly given by 

x{y) = n {4>,Y)(t)d(j), (9) 



and we will now calculate its statistical properties. 

Using equations (|3|), (|^) and (^), we see that the mean is given by 

= /^^(n'(0,f))p)^</># = B- {f{v))f, (10) 
JO 

where we have defined 

/(f) = n'(<A,)<A,2^(f ) (11) 

and where B denotes the average total brightness due to sources below our 
flux cut, i.e., 

B = [ n\(l))(t)d(l). (12) 



Since (A(r))j = (there is just as much positive as negative weak lensing 
from large-scale structure), the second term in equation (|lO|) vanishes, and 
we see that regardless of the flux cut, lensing has no impact on the average 
brightness of unremoved point sources. 

We now turn to the correlation function. Using (A(r))j = and equa- 
tions (^), (|^, (§) and (^), we obtain 

{x{v)x{v')) = j^^ j\{n'{ct>,v)n'[<P\v'))^)fckct>'dcPd<P' 

= ([5 -/(f)][i? -/(?')])/ + '5(f,?')^^^^ (13) 
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where 

CP'= / ' n'{(t))(t)'^d4> (14) 

is the familiar power spectrum of unresolved point sources in the absence of 
lensing that was derived in TE96. Using (/(?)) j = again, this reduces to 

(x(f )x(f')) = B^ + CP' + n'(0,)20,4(A(?)A(fO)/. (15) 

Since the first term is merely the familiar and uninteresting monopole, the 
new effect that we have computed corresponds to the third term. Thus the 
effect of lensing is to add to the original point source correlation function a 
new term which is the magnification correlation times n'(</>c)^</>c^. 



2.3 The effect on the point source power spectrum 

Let us expand the sky brightness x in spherical harmonics and investigate 
the statistical properties of the expansion coefficients 

a£m = J Y(*^{r)x{r)dn. (16) 

Assuming that the statistical properties of the magnification field A are 
isotropic (there is no need to assume that A is Gaussian), its corresponding 
expansion coefficients A^^ must satisfy 

i'^em^i'm') f = See'5mm'C^'^^ (17) 



rum. 



for some C™"^ which we will refer to as the magnification power spect 
Substituting equation (15) into 

(a|„a,,„,) = / / Y;^{r)Ye,^,{r){x{f)x{9'))dndfl', (18) 



we obtain 

{a}j^a£im') = See'Smm'Ce, (19) 
where the point source power spectrum is 

Ci = 47T6oe6ornB^ + + n'(</.,)2</,,4c77«5. (20) 

The first two terms are merely the monopole and the shot noise power (which 
is independent of £), derived in TE96, so the new effect due to lensing is 
given by the third term. 
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3 HOW IMPORTANT IS THE NEW EFFECT? 



3.1 Some useful approximations 

To clarify the relative importance of the last two terms in equation (|20|), it 
is convenient to express them in terms of quantities that are more directly 
linked to observations. For typical scenarios, the differential source count 
n'{4>) is a smooth function on a log-log-plot, which means that near the flux 
cut (pc, we can approximate it with a power law 

n'{<P)^n'i<p,)(Pi " (21) 



for some constant (3. To avoid the above-mentioned divergence of the total 
brightness, we must have a logarithmic slope /? < 2 at the faint end and 
/3 > 2 at the bright end. For instance, for radio sources at 1.5 GHz, the 
VLA FIRST survey gives the logarithmic slope /? ~ 1.6 at </)c ~ 1 mJy, 
steepening to /? ~ 2.5 at the bright end (White et al. 1996, TE96) just as 
one would expect if the brightest sources are at low redshifts z 1 where 
evolutionary and cosmological effects are negligible, and the same qualitative 
behavior is found at higher frequencies (Windhorst et al. 1985, 1993). Using 
equation (|2l|), we thus find that the expected number of sources above 
the flux cut is 

iV = 47r / n'{<P)d<P ^ —— </,,n'(0,) oc (t>c~^^-^\ (22) 

since the integral is dominated by sources just above the flux cut for which 
the power law fit is very accurate. Similarly, the integral in equation (p) is 
dominated by sources just below the cut, so using equations (|2l] ) and (2^) 
gives 

C"' - i ^] i ^) 4>c^ oc 0,3-/5. (23) 



3-/?/ V47r 



Finally, we can rewrite the lensing term in equation (20) as 



n'f 



^ {(5- 1)(3 - 0) c'^^acps 

oc ct^c-^^^-^^. (24) 
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Equation (|2^) showed that that power spectrum was a sum of three different 
terms. These are all plotted in Figure |2| as a function of the flux cut, using 
the source counts n'{(t)) from the VLA FIRST survey (White et al. 1996). 
The approximations above allow us to understand all qualitative features 
of this figure. Comparing equations ( p^ and (^), we notice that although 
the Poisson term C^'^ decreases if we remove more sources, the lensing term 
increases as long as /3 > 2, i.e., as long as we remove less than ~ 10^ sources 
from an all-sky survey, or ~ 4 sources per square degree. This corresponds 
to ~ 60 sources in a 4° x 4° field of the upcoming Very Small Array ( VSA) 
experiment, and to a 50 mJy flux cut at 1.5 GHz. 

3.2 The magnification power spectrum 

Equation (^) shows that when /3 ~ 2, the ratio between the lensing term 
and the conventional shot noise term C^^ is simply C™"^ times N/4tt, the 
number of removed sources per steradian. N/Att is typically 1, and we will 
now evaluate C^""^ . For any flat Friedman-Robertson- Walker cosmology, 
including non-linear clustering and an arbitrary source distribution n{xs), 
the magnification power spectrum is given by (Villumsen 1996) 

^mag ^ 'J2^^^2^-^ T'^ W^{xi){l + zfP{£/yd, z) dxi, (25) 

Jo 

where the function w is defined as 

w{xi) = I —n{xs) dxs- (26) 

Jxi Vs 

Here d = H^^c ^ 3000/i-^Mpc, P{k, z) is the conventional three-dimensional 
matter power spectrum at redshift z, and xd and yd are the comoving ra- 
dial and angular distances to the epoch z, respectively. Subscripts h, I and 
s refer to horizon, lens and source positions, respectively, and yis is the co- 
moving angular lens-source distance. These equations are valid on angular 
scales where the sky is approximately flat, i.e., ^ ^ 60, which is of course 
the case in our regime of interest. For the simple case where = 1, A = 0, 
the density fluctuations grow according to linear theory and all sources are 
located at some characteristic redshift Zg, this reduces to 

cr-=72A-/;(i-i-)V(A).., (.7) 

where P is the current power spectrum and Xg = 2 — 2/ \J1 + Zg. This is 
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Parameter 


Opt. 


Mid. 


Pess. 


Source redshift 


0.4 


1 


»1 


a 


0.2 





-0.3 




0.2 


0.7 


1 


Removed/sq deg 


0.5 


2 


4 




2.5 


2.2 


2 



Table 1: Parameters used in optimistic, middle-of-the-road and pessimistic 
scenarios. 



plotted in Figure ^ (dashed curves) for a standard CDM power spectrum 
(Bond & Efstathiou 1984) with h = 0.5, for the three cases described in 
Table 1. We label these scenarios "pessimistic", "middle-of-the-road" and 
"optimistic" , since they are intended to give an upper limit, realistic estimate 
and lower limit, respectively, as to how much of a problem the lensing effect 
will turn out to be for future CMB experiments. Apart from a factor of In 10, 
the magnification fluctuation on a given scale is essentially the area under 
a plotted curve out to that scale. For instance, the r.m.s. magnification 
fluctuations after smoothing on a scale of 5 arcminutes is about 15% for 
the pessimistic model, where the sources that have the dominant effect on 
the CMB are assumed to be at redshifts ^ 1. However, what matters for 
our purposes is of course not this r.m.s. magnification fluctuation, but C^""^ 
itself. In the pessimistic scenario, the maximum power is C^""^ ~ xl0~^, 
attained for £ ~ 10^, which ensures that the r.m.s. lensing fluctuations are 
below the standard Poisson fluctuations for any reasonable flux cut. 



3.3 The uncertain point source normalization 

This conclusion in independent of the overall normalization of the point 
source counts n'{4>), since Equation (24) shows that the ratio of our lens- 
ing effect to the shot noise power C^'^ is independent of this normalization. 
However, to assess the overall importance of both effects to future CMB 
experiments, accurate knowledge of the radio source counts at the frequen- 
cies where the observations are made is crucial. Unfortunately, the point 
source population between 20 and a few hundred GHz is still very poorly 
determined (see e.g. Gundersen et al. 1997), leaving open only the unattrac- 
tive option of extrapolating from lower frequencies. In Figure |^ we have 
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extrapolated the 1.5 GHz VLA FIRST data to 30 GHz by assuming that 
the temperature fluctuations (C^*)^/^(z^) oc z^~(^+"). To reflect the observa- 
tional uncertainty, we have allowed a to vary as shown in Table 1, with the 
middle-of-the-road estimate being that of TE96. 

4 CONCLUSIONS 

We have calculated the effect of weak gravitational lensing on the power 
spectrum of unresolved point sources. We found that this adds a new term, 
which is the magnification power spectrum times n' {(pc)'^ (pc'^ . The ratio 
between this term and the standard Poisson term is roughly (A^/47r)C™"^, 
independent of the (poorly known) point source normalization, so if we keep 
shrinking the point source region in Figure |^ by lowering the flux cut ad 
infinitum, the lensing term will eventually become the dominant problem. 
The lensing term initially increases as we remove more sources, reaching a 
maximum when ~ 4 sources per square degree are removed. 

Figure |l| shows that for high-resolution interferometric experiments op- 
erating at low frequencies, a fairly aggressive source removal scheme will be 
necessary, bringing the lensing effect close to this maximum strength. Effec- 
tive point source removal will be important for the future Satellite missions 
MAP and Planck as well. Fortunately, this strength does not exceed 12% of 
the standard Poisson fluctuations even in the worst-case scenario. 

In summary, although this effect should be taken into account in the 
detailed analysis of some future CMB experiments, it is likely to be a rel- 
atively small correction (like, e.g., the Rees-Sciama effect) and should not 
substantially reduce the accuracy with which cosmological parameters can 
be measured. 

The authors wish to thank Avi Loeb for suggesting this project and for 
helpful comments on the manuscript. This work was partially supported 
by European Union contract CHRX-CT93-0120, by Deutsche Forschungs- 
gemeinschaft grant SFB-375 and by NASA through a Hubble Fellowship, 
#HF-01084.01-96A, awarded by the Space Telescope Science Institute, which 
is operated by AURA, Inc. under NASA contract NAS5-26555. 
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Tegmark & Villumsen (1997) 




Figure 1: Where various foregrounds dominate. 
The shaded regions indicate where the various foregrounds cause fluctua- 
tions exceeding those of COBE-normahzed scale-invariant fluctuations, thus 
posing a substantial challenge to estimation of genuine CMB fluctuations. 
They correspond to dust (top), free-free emission (lower left), synchrotron 
radiation (lower left, vertically shaded) and point sources (lower and up- 
per right). The heavy dashed line shows the frequency where the total 
foreground contribution to each multipole is minimal. The boxes roughly 
indicate the range of multipoles i and frequencies v probed by various CMB 
experiments, as in TE96. 
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Figure 2: Dependence of radio source fluctuations on flux cut. 
The three terms that contribute to the point source power spectrum in 
equation ([20| ) are plotted as a function of the number of sources removed 
in an all-sky survey (only the number per steradian matters), based on the 
VLA FIRST source counts. These terms (solid curves) are the integrated 
brightness of the population, the edge effect at the flux cut and the shot noise 
fluctuations, respectively. The dashed curves show the lensing contribution 

1/2 

to for the pessimistic and middle-of-the-road scenarios, at the worst- 
case values. 
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Figure 3: The power spectra are plotted for the CMB (heavy Une), the 
lensing magnification field (the three top lines), discrete point sources at 30 
GHz (the three straight lines) and the lensing effect on them (the two bottom 
lines), all for a standard CDM model. The solid, dashed and dotted lines 
correspond to the pessimistic, middle-of-the-road and optimistic scenarios, 
respectively. 
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